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1. Introduction 

We live in a time when Standard Model (SM) has been vindicated in a number of experiments 
and at the same time has left us with a feeling that it is incomplete, which makes the search for 
physics beyond SM a holy grail for most of the particle physics community. We have absolutely 
no idea what the manifestation of the new physics will be; this forces us to consider all sorts of 
theoretically motivated new physics options. Since one does not have any experimental signal that 
definitely points to new physics, the best one can do is to constrain the parameter space of the new 
physics models. These constraints come mostly from experimental data (including astrophysical 
ones) but sometimes from theoretical considerations too []. 

As long as one does not produce the new particles directly, one has to look for their indirect 
effects in low-energy observables. All low-energy data are more or less consistent with the SM, 
taking into account the experimental and theoretical errors and uncertainties. Thus, one can look 
at those observables which may be explained by SM; the trick is to maximize the width of the 
window for new physics and put constraints on the parameter space of such models. Alternatively, 
one can look at those observables which are absolutely forbidden or highly suppressed in SM so 
that one does not expect any signal; here even one event will signal new physics and in the absence 
of any event, the parameter space for new physics may be constrained from the experimental upper 
bounds. 

Among the new physics options that people consider, supersymmetry (SUSY), with all its 
variants, is the most popular one M, In the minimal and some non-minimal versions of SUSY, 
the action is so taken as to conserve the R quantum number defined as R = (— \^3B+L+2S ^ where 
B, L and S stand for baryon number, lepton number and spin of the field respectively. This 
ensures R = +1 for all particles and R = — 1 for all superparticles, and conservation of R imply 
that superparticles must occur in pair in all allowed Feynman vertices. However, R-parity is a 
discrete symmetry imposed by hand (to make the parameter space of the model more restricted 
and tractable) and one can write a R-parity violating superpotential of the form 

W = \\ l]k UL 3 El + X' ijk LiQjDl + ±\1 jk UfDp c k (1) 

which does not violate any gauge symmetry. The factors of 1/2 take care of the fact that A and A" 
couplings are antisymmetric in its first two indices. Such R-parity violating terms can be motivated 
from some of the grand unified theories Q. In eq. ([[]) L, Q, U, D and E denote, respectively, 
SU(2)l doublet lepton and quark superfields, and SU{2)l singlet up, down and charged lepton 
superfields, and k are generation indices. Of course, B and L are both violated, and to forbid 
proton decay, one has to keep either L-violating A and A' terms or i?-violating A" terms, but 
not both. There is also a bilinear R-parity violating term of the form e^Lj^j which has a lot 
of interesting phenomenology ||, including possible leptonic flavor violation H, but we will not 
consider that term in the present paper. 

With 45 new couplings (and all of them can theoretically be complex) the phenomenology is 
immensely richer, but at the same time less predictive. There is, however, one major point to 
be noted: R-conserving SUSY can affect low-energy observables through loop effects and hence 
can hardly compete with SM effects (except in some of the cases where the SM process itself 
is loop-induced or, even better, forbidden); R-parity violating (RPV) SUSY, on the other hand, 
can show up in tree-level slepton or squark mediated processes which can successfully compete 

1 For example, from the consideration of the stability of the scalar potential. 
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with the SM. This also ensures that for comparable coupling strength, RPV SUSY amplitudes are 
generally orders of magnitude higher than the R-conserving SUSY amplitudes. 

The individual RPV couplings have been constrained from various low-energy processes ||, 
and upper limits on some of the product couplings, including their phases, have also been found 
|], [7|, ||, ^, [l0|, O]. Often one finds that the product coupling is much more constrained than the 
direct product of the upper bounds of the individual couplings f| 

In this paper we find the upper limits on the products of some of the L-violating A and A' 
type couplings coming from rare r, Bj and B s decays. The leptonic flavour violating processes are 
forbidden in the SM, and the only contributing amplitude comes from RPV SUSY. The expected 
branching ratios (BR) of leptonic flavour conserving AB = 1 processes within SM are so much 
below the experimental numbers (except B — ► K^£ + £~) that one can safely ignore the SM effects, 
as well as the R-conserving SUSY effects, to put bounds on the RPV couplings. Most of these 
decay modes are also theoretically clean and free from any hadronic uncertainties which plague 
nonleptonic decays. (The exception is again the semileptonic process 

There are three major sources of uncertainty, however. First, the decay constants of the 
neutral mesons, particularly that of 77, B^ and B s , are yet to be cleanly determined. Fortunately, 
the bounds where these decay constants are relevant just scale with their values; this will be 
discussed in the appropriate section. Second, the B — > K/ir formfactors allow a slight theoretical 
uncertainty. We use the BSW formfactors with oscillator parameter oj = 0.5 GeV |l2|]. Lastly, the 
current masses of the light quarks cause the maximum uncertainty in AA' type couplings; again, 
the bounds scale with the quark masses. The individual bounds in most of the cases are fairly 
weak and the bounds on the product couplings that we find are sometimes one to two orders of 
magnitude better than the existing ones. Some of these product couplings have been considered 
earlier |@, ||; these we update with better experimental numbers. These updated numbers, too, 
are vastly improved. We have not discussed leptonic and semileptonic D° decays since the bounds 
are much weaker than those which one gets from — ► e conversion ||. This in turn implies that 
such decay signals of D°, if observed in present and upcoming colliders, imply some new physics 
but not RPV SUSY. 

When we say that these bounds are robust, we of course not only mean to alleviate the 
theoretical uncertainties of the SM. The fact that leptonic flavour-changing processes are absolutely 
forbidden in the SM imply that these bounds stand no matter what the phases of these RPV 
couplings may be. The same is true for all purely leptonic B decays, but not if there are competing 
SM amplitudes (as in B -> Kt + l~), even if one takes the lowest possible SM number and saturates 
the experimental data by RPV contribution, which is the standard practice. The reason is simple: 
the two amplitudes are coherent and the bounds depend on the phase of the RPV couplings, while 
the standard prescription is true for incoherent amplitudes only. This we show explicitly for the 
B — > K£ + £~ decays. One should note here that the most conservative bounds, not necessarily the 
true ones, come from such incoherent amplitude summation. 

There is another point that we like to emphasize. Though it is true that observation of the SM- 
forbidden decays would be a definite signal of new physics, to show convincingly that it is RPV 
SUSY one needs to find some correlated signals (e.g., enhancement of BRs) in different decay 
channels. We briefly discuss how this can be done for explicit RPV models (spontaneous RPV 



models have been discussed in [13|, including bounds coming from mesonic and leptonic flavor 



violating decays). If one finds such correlated signals, it will be an almost definite RPV signal, 

2 Constraints coming from ArriK, Ams, ft — > 3e, fj, —* e conversion, etc. fall in this category. 
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without the direct observation of superparticles. Such signals occur from the fact that same four- 
Fermi operators lead to different final states. Among leptonic and semileptonic B decays, one may 
mention the correlated channels (i) Bd — ► e^i and Bd — > Trefi; (ii) B s — ► if£j and Bd — ► KifiJ. 
More correlated channels are to be found in nonleptonic B decays, which will be discussed in 



a future paper [14]. In the absence of any correlated decay channel, the next best thing is to 
observe the decay distribution of the final state particles, since RPV SUSY has a different Lorentz 
structure from that of the SM. One may, for example, study the angular distribution of final state 
leptons in B — > Kl + l~ decays — the crucial fact is that tree-level RPV has a Lorentz structure of 
the form (V — A) (g) (V + A) . We do not go into any detailed discussion of this issue in the present 
work. 

We consider only nonzero L- violating A and A' type couplings. Though it is true that leptonic 
r decays can be mediated by A type couplings alone, the individual bounds on these couplings 
are much tighter than one may hope to get from such decays. Thus, r — > 3 lepton (there can be 
six different combinations) processes give a fairly weak bound (~ 0(1)) on the relevant A- type 
couplings. We do not consider such processes further; if they are observed in near future, their 
explanation must lie somewhere else. It is obvious that A" type couplings cannot mediate leptonic 
and semileptonic decays. 

To construct four-Fermi operators from A and A' type couplings that mediate such semileptonic 
and leptonic r and B decays, one needs to integrate out the squark or the slepton propagator. 
Both the couplings coming in the product may be complex; one is free to absorb the phase of one 
coupling in the sfermion propagator but the other remains, making the overall coupling responsible 
for the process a complex one in general. However, since all these processes are one-amplitude 
(there is no SM counterpart) no scope of CP-violation exists; one only observes the nonzero 
branching ratios. By the same argument, we can take all couplings to be real without any loss 
of generality. The only exception to this statement, viz. B — > K(*)£ + £~ , will be dealt in proper 
place. 

The paper is arranged as follows. In the next section, we discuss the formalism for, and the 
bounds coming from, semileptonic r decays. Section 3 discusses leptonic Bd and B s decays, while 
section 4 is on semileptonic decays of these mesons. We summarize and conclude in section 5. 



2. r decays 

All the processes that we consider involve terms in the RPV Hamiltonian with two leptons 
and two quarks as external fields. The Hamiltonian can be written as 

H 9 = -A jk i m (£j{l + 7 5 )4)(^(1 - 7sM) + 2 5 ife/m(^7 M (l - 7s)^)(&m(! + 7s)4) 

-\Cjklm{Tj^{l - 7 5 )^)(^7m( 1 + 75)^m) + h.c. (2) 

where 

3 A* A' 3 A'* A' 3 A'* A' 

a _ \ " A ijk A ilm _ R _ \ ~> A jik A lim _ p, _ \ ^ A jki A lmi /r,\ 

t— 4m,->2 r— ' 4m, r 2 ^— f 4m ,- 2 

8=1 v i 1=1 u L t 1=1 d Ri 

To put bounds, we will assume only one of the A, B or C terms to be nonzero so that there 
is no interference effect between different RPV couplings. We will also take only one sfermion 
generation index i to be nonzero at a time. 
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Note that if the final state consists of two down-type quarks, both A and B terms may 
contribute, whereas for two up- type quarks, only the C term comes. For mesons like 7r°, p° or r\ 
in the final state, all the three terms may be important. 

One thing that we do not consider is the running of the RPV couplings between the sfermion 
scale and the low-energy scale. The corrections are electroweak in origin and can be safely ne- 
glected. The only QCD correction may occur between the two quark fields; for r decays, this pair 
hadronizes, absorbing all such uncertainties in the decay constant. The same is true for leptonic 
B or D decays. For semileptonic meson decays the formfactors are supposed to take care of these 
short-distance corrections. However, this effect is important when we have four quarks as external 
fields and may contribute a multiplicative factor of ~ 2 to the effective Hamiltonian . 

The generic process r — > i + M is lepton-flavour violating and does not occur in the SM. 
Strong experimental upper limits exist on at least fourteen modes that we consider here: £ = e,p, 
and M = tt, p, ry, K°,K° ,K *, (/). Such modes are fairly clean from a theoretical point of view. 
In RPV SUSY, all such processes can occur with squark or sneutrino propagator mediating the 
decay. Since the squark is very heavy, we do not consider any QCD effect that may take place 
between the squark and the final state quarks. 

The only uncertainty, albeit small, appears in the decay constants of the neutral mesons. Our 
values for the decay constants are (in GeV) fl6|] : 

U = 0.132, f p = 0.216, f K = 0.161, f K . = 0.214, / = 0.237. (4) 

The decay constants for rj (and rf) are obtained from the decay constants of the octet and singlet 
mesons fg = 1.34/^ and /i = 1.10/„- by a rotation: 

fv = fi = h cos 0/V6-/i sin 0/V3, 

/• = -2/ 8 cos e/VE-fi sin e/Vs, 

# = / 8 sin0/\/6 + /ico S 0/\/3, 

/*, = -2/ 8 sin e/Ve + /i cos e/Vs, (5) 

where the angle 9 is estimated to be about —22°. The current quark masses are taken to be 
m d = 10 MeV and m s = 200 MeV. 

Following the standard practice, we assume only one product of RPV coupling to be nonzero at 
a time. This eliminates the need to consider their phases and signs; without any loss of generality, 
we can assume all products to be real and positive since only the square of the absolute magnitude 
of the product enters in the expression for the decay width. If the final state quarks are of charge 
—1/3, the mediating squark must be a 'left-handed' up-type one, and if the quarks of charge +2/3, 
the mediating squark is a 'right-handed' down-type one. Note that if two or more products are 
simultaneously nonzero, there can in principle be an interference effect for M = tt, p and 77, and 
it becomes imperative to consider the signs and the phases of the product couplings. We neglect 
this complexity. 

The experimental 90% CL upper limits on the BRs of various r decay modes are taken from 

Br(r -> evr) < 3.7 x 10" 6 Br(r -» /ivr) < 4.0 x 10" 6 

Br(r -► erf) < 8.2 x 10~ 6 Br(r -► prf) < 9.6 x 10~ 6 

Br(r -> eK°) < 1.3 x 10~ 3 Br(r -» pK°) < 1.0 x 10~ 3 
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Br{r -> ep) < 2.0 x 10~ 6 
Br(r -> e#) < 6.9 x 10~ 6 
Br(r -» eK*°) < 5.1 x 10~ 6 
5r(r -> eK*°) < 7.4 x 10~ 6 



Br(r pp) < 6.3 x 10" 
Br(j -> p4>) < 7.0 x 10" 



Br(r 
Br(r 



pK 



*(K 



pK 



*0\ 



< 7.5 x 10" 

< 7.5 x 10" 



(6) 



The nonzero A A' type couplings can mediate only r — > ^ + P type decays where P is a generic 
pseudoscalar meson. That the production of vector mesons is forbidden is evident from the Lorentz 
structure of the corresponding four-fermi hamiltonian. The decay width can be written as 



where 



T(T^£i + P[=q-q k }) 



C(mi,m 2 ,m 3 ) 



(m 2 + m|. — mp)C(m T , m£.,mp)Fp 



128vrm3m 4 



njk I 



mf + 771-2 + m 3 



2mfm3 — 27712^3, 



Pp 



4/? 



(m 9i + m 



(7) 

(8) 
(9) 



and 77i denote the slepton mass. Note that due to our definition of the decay constants, the above 
expression is to be multiplied by a factor of 1/2 only if there is a ir° in the final state. The 
combination A^A^ can be replaced by A n 3jA^ fcj - which will generate the same decay and hence 
the same bound applies to both these combinations. 
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Table 1: Bounds on AA' type products from r — > t + M decays. 



The bounds are listed in table 1. Note that most of the numbers are at the same order of 
magnitude as the previous bounds coming from direct product of the individual couplings, though 
some have been improved^. We show only the best numbers; for example, the decay r — > e + ir° 
puts a weaker bound on A23iA 211 than that coming from r — > e + p, and hence is not shown 
separately. Here, and in all other cases, we take all squarks and sleptons to be degenerate at 

3 Previous bounds indicate the numbers coming from direct product of individual bounds, or, in certain cases, 
bounds on the product coming from a different process. The bounds which are updated have not been considered 
as previous bounds. 
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100 GeV, and the bounds scale in a simple way: (m-^/100 GeV) 2 . For squarks of the first two 
generation this number is not allowed but used just as a benchmark value; for the lighter stop 100 
GeV is still allowed, and light sbottom is of current phenomenological interest. Thus, as far as 
squark-mediated processes are concerned, the more realistic bounds for the first two generation of 
squarks should be the number quoted in our tables multiplied by a factor of ~ 10. 

Though we have relied on published experimental numbers only to obtain the bounds, let us 
also quote the Belle numbers for the modes r — > eK° and r — ► fiK° 

Br(r -» eK°) < 1.8 x 10~ 6 , Br(r -> pK ) < 1.8 x 10" 6 . (10) 

If we use the Belle data, all the eight AA' type products in table 1 where the final state is either eK° 
or fiK° have upper bounds of 1.7 x 10 -3 , which is better or compatible to the previous numbers. 

The generic couplings B and C (eq. ||) mediate the decay r — > £ + M where M can be either 
pseudoscalar or vector. The expression for the decay widths in these two CcLSGS £1X6 ELS follows: 

' C(m T , m 

512vrm 4 m3 \*3nk*inj\ 



9 j.pr-7i-«n / p g(m T , m ti , m P )P (m T , m ti , m P ) , , 2 
1 (r — > ti + f[= qjqkl) = — — — — 3 |A 3n . fc Ai„.,-| 



r(r -» £j + K[= qjqkl) - 5127rm 4 m3 \ X 3nk\nj\ (H) 

where C(m T ,m£.,my) is defined in eq. (||) and 

P (*,y,2) = (x 2 -y 2 ) 2 -z 2 (^ 2 + y 2 ) 

y (x,y,z) = zV + y 2 -z 2 ) + (* 2 -y 2 ) 2 -z 4 (12) 

The combination Ag nfc A^ n:; - appears if both qj and qt are down-type quarks. In this case the 
mediating squark is u n L- If t ne final-state quarks are up- type, the combination that appears 
is X' 3 j n X' ikn , with rest of the formula remaining unchanged. Again, the expressions are to be 
multiplied by a factor of 1/2 if one has 7r° or p° in the final state. 

We list only the best bounds coming from these processes in table 2. Note that the best 
bounds always come from those decays where a vector meson is involved (if we do not consider 
the Belle data in eq. [U]) in the final state; this is solely due to the better limits on the decay 
modes. Since some of the A' type couplings have weak individual bounds, most of our bounds on 
the product couplings are one to two orders of magnitude improvement over the previous bounds. 

Consideration of Belle data in eq. (10) shows that the couplings A'^A^ and A^A^ both 
have upper bounds of 2.3 x 10~ 3 coming from r — > eK° and r — > fj,K° respectively, all other 
bounds remaining unchanged. 

Let us, at this point, highlight certain features of the analysis: 

• Mere observation of a single event in any of the decay modes will signal new physics. 

• However, in a dedicated tau-charm factory, one may hope to observe more events in different 
channels if the RPV couplings are close to their present bounds obtained in tables 1 and 2. 

• With only one nonzero A'A' product, one should observe signals in different modes corre- 
sponding to the same quark-level subprocess: e.g., r — * eTT,erj,ep all should show some 
anomalous behaviour. This is another example of correlated channels and is important to 
establish the nature of the new physics. If no such signals are observed in the r — > p + M 
channels, flavour-specific nature of the new physics will all the more be established. 
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Table 2: Bounds on A'A' products from r — > 1 + M processes. Only the best bounds are shown. 



• If signals are observed in the pseudoscalar channels but not in the vector channels, that 
will probably indicate the presence of a AA' type coupling. If the opposite happens, RPV 
explanation of new physics will be more difficult to sustain; probably one has to invoke 
cancellation between different RPV contributions. 

• Finally, note that though the individual A or A' type couplings are L-violating, the overall 
four-fermi hamiltonian conserves L (it violates leptonic flavour). Thus, one cannot explain 
processes like r~ -> t r M^M 2 with RPV. 



3. Leptonic Bd and B s decays 

The leptonic flavour-violating decays B^ jS — ► ^t^-J {i 7^ j) are forbidden in the SM, and 
flavour-conserving decays (i = j) are so suppressed (except for t = r which we do not consider 
anyway) that we can take them to be almost forbidden to a very good extent. Thus, the entire 
amplitude, if nonzero, is solely due to new physics. In RPV models, both slepton-mediated AA' 
type and squark-mediated A'A' type interactions can cause such purely leptonic decays. As already 
stressed, the bounds are robust in the sense that they are free from any theoretical uncertainties 
(except for the decay constants of Bd and B s ), and do not depend on the phase of the RPV 
couplings. 

The decay width of Bd tS — >■ is given by 



/ 2 



nB qi - = lte ^3 C ( M B gi , m tl , m tm )Px (M Bqi , m k ,m e J\ \ nlm X' ni3 \ 2 (13) 

Bo- 



ot 



f 2 

T{B qi qi-J = *l M3 C(M Bqi , m ei ,m e JP 2 (M Bqi , m tr m im )\X' lni X' mn3 \ 2 . (14) 
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Here Pi and P2 are given by 

n / \ 4 22 22 

Fi{x,y,z) = x — x y — x z , 

P 2 (x,y,z) = x 2 (y 2 + ^ 2 )-(y 2 -z 2 ) 2 (15) 

and the C-function is defined in eq. (|8|). The generic slepton or squark mass is denoted by to. 
Note that the expression A ra ; m A^ i3 can be replaced by A ram ;A^ 3 j in eq. (|i"3"l) and A; ni A^ 7in3 can be 
replaced by A^gA^j in eq. fll4|). Since we consider only one combination to be present at a time, 
the same bound applies to all such combinations. 



AA' 


Process 


Bound 

(xlO 5 ) 


Previous 
bound 


AA' 


Process 


Bound 

(xlO 5 ) 


Previous 
bound 


(121) 


(113) 


B d - 


-> fjpe* 


2.3 


1.0 x 


10" 


-3 


(121) 


(123) 


B s - 


-> fjPe* 


4.7 


2.1 x 10" 


-3 


(121) 


(131) 


Bd - 


-> H*e^ 


2.3 


9.3 x 


10" 


-4 


(121) 


(132) 


Be - 


-> ) u ± e =F 


4.7 


1.4 x 10" 


-2 


(121) 


(213) 


Bri - 


-* e + e~ 


1.7 


2.9 x 


10" 


-3 


(121) 


(223) 


B s - 


-> e + e — 


14 


1.0 x 10" 


-3 


(121) 


(231) 


B d - 


-> e + e _ 


1.7 


8.8 x 


10" 


-3 


(121) 


(232) 


B s - 


-> e + e — 


14 


7.3 x 10" 


-3 


(122) 


(113) 


Bh - 

(I 


•* u u~ 


1.5 


5.4 x 


10" 


-3 


(122) 


(123) 


B« - 


-> u u~ 


2.7 


2.1 x 10" 


-3 


(122) 


(131) 


B d - 


-> fl + 'fjT 


1.5 


9.3 x 


10" 


-4 


(122) 


(132) 


B s - 




2.7 


1.4 x 10" 


-2 


(122) 


(213) 


B d ~ 


-> ^e* 


2.3 


2.9 x 


10" 


-3 


(122) 


(223) 


B s - 




4.7 


1.0 x 10" 


-2 


(122) 


(231) 


B d ~ 


^^e T 


2.3 


8.8 x 


10" 


-3 


(122) 


(232) 


B s - 




4.7 


2.7 x 10" 


-2 


(123) 


(113) 


B d ~ 


+ ^t^ 


62 


1.0 x 


10" 


-3 


(123) 


(131) 


B d ~ 




62 


9.3 x 10" 


-4 


(123) 


(213) 


B d ~ 


-» e ± r^ 


49 


2.9 x 


10" 


-3 


(123) 


(231) 


B d ~ 


-> e ± r^ 


49 


8.8 x 10" 


-3 


(131) 


(113) 


B d ~ 


-» e ± r^ 


49 


1.3 x 


10" 


-3 


(131) 


(131) 


B d ~ 


-> e ± r^ 


49 


1.1 x 10" 


-3 


(131) 


(313) 


B d ~ 


-» e + e~ 


1.7 


6.8 x 


10" 


-3 


(131) 


(323) 


B s - 


-» e + e~ 


14 


3.2 x 10" 


-3 


(131) 


(331) 


B d ~ 


-> e + e~ 


1.7 


2.8 x 


10" 


-2 


(131) 


(332) 


B s - 


-> e + e~ 


14 


2.8 x 10" 


-2 


(132) 


(113) 


B d ~ 




62 


1.3 x 


10" 


-3 


(132) 


(131) 


B d ~ 




62 


1.2 x 10" 


-3 


(132) 


(313) 


B d ~ 




2.3 


6.8 x 


10" 


-3 


(132) 


(323) 


B s - 




4.7 


4.3 x 10" 


-3 


(132) 


(331) 


B d ~ 




2.3 


2.8 x 


10" 


-2 


(132) 


(332) 


B s - 




4.7 


2.8 x 10" 


-2 


(133) 


(313) 


B d ~ 


-►eVF 


49 


3.6 x 


10" 


-5 


(133) 


(331) 


B d ~ 


^eV 


49 


2.7 x 10" 


-3 


(231) 


(213) 


B d ~ 


-►eVF 


49 


4.1 x 


10" 


-3 


(231) 


(231) 


B d ~ 


^eV 


49 


1.3 x 10" 


-3 


(231) 


(313) 


B d ~ 




2.3 


4.2 x 


10" 


-4 


(231) 


(323) 


B s - 


^^e T 


4.7 


3.6 x 10- 


-2 


(231) 


(331) 


B d ~ 




2.3 


2.2 x 


10" 


-2 


(231) 


(332) 


B s - 




4.7 


3.2 x 10- 


-2 


(232) 


(213) 


B d ~ 




62 


4.1 x 


10" 


-3 


(232) 


(231) 


B d - 




62 


1.3 x 10" 


-2 


(232) 


(313) 


B d ~ 




1.5 


1.3 x 


10" 


-2 


(232) 


(323) 


B s - 




2.7 


3.6 x 10- 


-2 


(232) 


(331) 


B d ~ 




1.5 


3.5 x 


10" 


-3 


(232) 


(332) 


B s - 


+ ia + (a- 


2.7 


3.2 x 10- 


-2 


(233) 


(313) 


B d ~ 




62 


7.7 x 


10" 


-3 


(233) 


(331) 


B d ~ 




62 


3.2 x 10" 


-2 



Table 3: AA' bounds from leptonic Bd and B s decays. 



We present our numbers for AA' couplings in table 3 and A' A' couplings in table 4. Our input 

'0 CL upper bounds on the BRs of the following modes [17, 



parameters are the experimental 



Br(B d - 


-» e+O < 8.3 x 10" 


-7. 


Br(B d - 


-> /x + /i") < 6.1 x 10~ 7 ; 




Br(B d - 


-» e ± /x T ) < 1.5 x 10" 


-6. 


Br(B d - 


-> e ± r T ) < 5.3 x 10" 4 ; 




Br(B d - 


+ p^r?) < 8.3 x 10" 


-4. 


Br(B s - 


■> e + e~) < 5.4 x 10~ 5 ; 




Br(B s - 


* Ai + /x") < 2.0 x 10" 


-6. 


Br(B s - 


-> e ± Ai =F ) < 6.1 x 10" 6 . 


(16) 



Furthermore, we take the decay constants of both Bd and B s to be 200 MeV; the bounds just 
scale as (/ Bds /200 MeV) 2 . 
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We also have Belle numbers for some of these modes 



[20] 



Br(B d e + e") < 6.3 x 10" 7 ; Br{B d fi + ^) < 2.8 x 10~ 7 ; Br{B d e ± /x =F ) < 9.4 x 1(T 7 . 

(17) 

If one takes these bounds, which are yet to be published, into account, some of the numbers in 
tables 3 and 4 get modified. In table 3, the AA' combinations coming from B d — ► e + e~ have upper 
bounds of 1.5 x 10~ 5 instead of 1.7 x 10~ 5 ; those coming from B d — ► /U + /U _ and .B^ — > /z^e^ 
have upper bounds of 1.0 x 10~ 5 and 1.8 x 10 -5 respectively. In table 4, the bounds on the A'A' 
combinations coming from B d — ► /U + /i~ and — > /x^e^ are modified to 1.4 x 10~ 3 and 3.7 x 10 -3 
respectively. Obviously, these bounds scale as the square root of the upper bounds on respective 
branching fractions. 



A'A' 


Process 


Bound 


Previous 


A'A' 


Process 


Bound 


Previous 










(xlO 3 ) 


bound 










(xlO 3 ) 


bound 


(111) 


(213) 


B d - 




4.7 


3 x 10" 5 


(111) 


(313) 


B d - 




5.9 


5.7 x 10~ 5 


(112) 


(213) 


B s - 




9.6 


1.2 x 10" 3 


(113) 


(211) 


B d - 


-> /i ± e T 


4.7 


1.2 x 10~ 3 


(113) 


(212) 


B s - 




9.6 


1.2 x 10" 3 


(113) 


(311) 


B d - 


-> e±T* 


5.9 


2.3 x 10" 3 


(121) 


(223) 


B d - 




4.7 


9.0 x 10~ 3 


(121) 


(323) 


B d - 


^eV 


5.9 


22.4 x 10~ 3 


(122) 


(223) 


B s - 




9.6 


9.0 x 10~ 3 


(123) 


(221) 


B d - 


-> f^e* 


4.7 


7.7 x 10" 3 


(123) 


(222) 


B s - 




9.6 


9.0 x 10~ 3 


(123) 


(321) 


B d - 




5.9 


22.4 x 10" 3 


(131) 


(233) 


B d - 




4.7 


2.9 x 10~ 3 


(131) 


(333) 


B d - 


^eV 


5.9 


7.6 x 10~ 3 


(132) 


(233) 


B s - 




9.6 


42 x 10~ 3 


(133) 


(231) 


B d - 


-> l^e* 


4.7 


0.3 x 10~ 3 


(133) 


(232) 


B s - 


-> n^e* 


9.6 


0.8 x 10~ 3 


(133) 


(331) 


B d - 


^eV 


5.9 


0.6 x 10~ 3 


(211) 


(213) 


B d - 


-> fi + fi- 


2.1 


3.5 x 10~ 3 


(211) 


(313) 


B d - 


-» /x±tT 


7.3 


6.5 x 10~ 3 


(212) 


(213) 


B s - 


-> h+ij,- 


3.9 


3.5 x 10~ 3 


(213) 


(311) 


B d - 


-» ^r^ 


7.3 


6.5 x 10~ 3 


(221) 


(223) 


B d - 




2.1 


1.4 x 10~ 3 


(221) 


(323) 


B d - 


/i ± r^ 


7.3 


93.6 x 10" 3 


(222) 


(223) 


B s - 




3.9 


2.7 x 10~ 3 


(223) 


(321) 


B d - 


+ /x±tT 


7.3 


0.109 


(231) 


(233) 


B d - 




2.1 


27 x 10~ 3 


(231) 


(333) 


B d - 


-» ^r^ 


7.3 


81 x 10~ 3 


(232) 


(233) 


B s - 




3.9 


84 x 10" 3 


(233) 


(331) 


B d - 




7.3 


67.5 x 10" 3 



Table 4: A'A' bounds from leptonic B d and B s decays. 



The present e + e~ B factories should improve these bounds by one order of magnitude at the 
end of their run, if such modes are not observed. The hadronic machines or high-luminosity e + e~ 
machines like the projected SuperBaBar [21| should see a large number of such leptonic decays if 
the actual values of the couplings are anywhere near the present bounds. Since the same couplings 
cause both B d — > and B — ► nlilj (the same is true for B s — ► elkij and B — ► Klitj) decays, a 
simultaneous signal is expected. 



4. Semileptonic decays: B — ► Ke+e ,B — > Kfi+n ,B — > Ke^fj,^, B — ► 7re 



All three collaborations CLEO, BaBar and Belle have set upper limits on the BRs of the 
abovementioned semileptonic modes (and also modes with a vector meson in the final state) 

Ke + e~ at less than 3a) have been 



||, H g§ |25|]. In fact ' the modes B -► Kfi + n~ (and B 
observed by Belle [24]. The present status is as follows: 



Br(B Kl + T) < 0.6 x 10~ 6 (BaBar), 1.49 x 10~ 6 (CLEO) (t = e/y) 
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Br(B - 
Br(B - 
Br(B 



> Ke + e~) 
-> he ± p T ) 



+0.32- 
0.24- 



0.09 n, 
0.11 J 



(0.481 

(o.99t^r i:5 



x 10~ 6 (Belle) 
x 10~ 6 (Belle) 



< 1.6 x 10" 



0.14 

(CLEO), 



(h = K/tt) 



(18) 



These numbers are at the same ballpark as the SM expectations. Unfortunately, the SM 
expectations are not precise; at least three different groups quoted three different numbers, which 
are all mutually compatible, but differ in their upper and lower limi ts p6| , p7| , pEj ]. For example, 
the predicted BR for the mode B -» Kp + p~ is (i) (0.57±g;ig) x 10~ 6 |f; (ii) (0.33 ±0.07) x lO" 6 
|27|; (iii) (0.42 ± 0.09) x 10~ 6 @. As we know, the bounds on new physics depend sensitively 
on the theoretical predictions. Moreover, the situation is complicated since we should not neglect 
the SM amplitudes (apart from the lepton flavour- violating cases). 

We take a compromising approach: in table 5, we quote the bounds on the relevant RPV 
couplings assuming (i) that the SM expectation is at its lowest possible value predicted by pq ]; (ii) 
that the experimental number is at its highest possible value; (iii) that the difference is saturated 
by RPV contribution; (iv) and, most important of all, that the total SM amplitude and the RPV 
amplitude add incoherently. This is certainly not true; for that, in table 6, we show how the 
bounds relevant for a particular decay mode (viz., B — > Kp + fi~) change if (i) the RPV amplitude 
adds constructively to the SM one; (ii) the RPV amplitude adds destructively with the SM one (for 
this, we take the highest possible SM prediction and the lowest possible experimental number); 
and (iii) if we use the theoretical numbers in [^7]] or [28]. We draw the attention of the reader 
to the huge fluctuation of the RPV bounds. It is, however, heartening to note that the most 
conservative bounds come from incoherent amplitude summation, so that the bounds quoted in 
table 5 are robust as far as the phases in the RPV couplings are concerned. This analysis also 
shows that to get any meaningful signal of new physics from these modes, which is certainly 
feasible in not-too-distant future, one must minimize the theoretical as well as the experimental 
uncertainties. Another signal of new physics, which we do not investigate here, is the forward- 
backward asymmetry of the final state leptons. 



A'A' 


Final 

state 


Bound 


A'A' 


Final 

state 


Bound 


(111)(213) 




5.0 x 


10" 


-4 


(113)(211) 


n^e^ 


5.0 x 


10" 


-4 


(121)(223) 


irfi^e^ 


5.0 x 


10" 


-4 


(123)(221) 


TTjJL^e^ 


5.0 x 


io- 


-4 


(131)(233) 


irfi^e^ 


5.0 x 


io- 


-4 


(133)(231) 


TTjJL^e^ 


5.0 x 


io- 


-4 


(112)(213) 


K^eT 


3.3 x 


io- 


-4 


(113)(212) 


K^eT 


3.3 x 


io- 


-4 


(122)(223) 




3.3 x 


io- 


-4 


(123)(222) 




3.3 x 


io- 


-4 


(132)(233) 


K^e^ 


3.3 x 


io- 


-4 


(133)(232) 


K^e^ 


3.3 x 


io- 


-4 


(212)(213) 


Kp+p- 


9.7 x 


io- 


-5 


(222)(223) 


Kp+p- 


9.7 x 


io- 


-5 


(232) (233) 


Kfi+fT 


9.7 x 


io- 


-5 


(112)(113) 


Ke + e- 


9.6 x 


io- 


-5 


(122)(123) 


Ke + e~ 


9.6 x 


10" 


-5 


(132)(133) 


Ke + e~ 


9.6 x 


10" 


-5 



Table 5: A'A' bounds from semileptonic Bj decays. 



For our analysis, we take the most stringent experimental numbers, viz., the numbers quoted 
by BaBar for B — > K£ + £~ (£ = e, n) and the numbers quoted by CLEO for B — > Kep, irep. The 
SM prediction for the first two modes can be as low as 0.47 x 10~ 6 [2C] and zero for the last two; 
with incoherent addition of SM and RPV amplitudes, one gets the bounds on A'A' type products 
as shown in table 5. We do not show the AA' bounds as they are at least one order of magnitude 
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weaker than those obtained from leptonic and B s decays, apart from being much less robust. 
We have also used the BSW formfactors Jl^] for B — > K and B —* ir transitions and taken Fq = F\ 
with minimal loss of accuracy: F^ -+7r (0) = 0.39 GeV, F^ K (0) = 0.42 GeV. 

Though the bounds in table 5 appear to be better than those in table 4, we warn the reader 
that they depend sensitively to the experimental number as well as the scheme for theoretical 
prediction. This is evident from table 6. Thus, any future analysis with these product couplings 
should use the numbers quoted in table 4. 

Is it possible to find other channels where one finds such RPV signals? If only one product 
coupling is nonzero, there are very few product couplings constrained here that may contribute 
to nonleptonic B decays. Such couplings are: (i) A 211 A' 2 i3, mediating b — > uud and b — > ddd, and 
hence B — > tttt (this coupling also affects the B^-B^ box diagram and thus the mixing-induced CP 
asymmetry coming from decays); (ii) A 212 A 21 3 and A' 112 A' 113 , mediating b — > uus and b — ► dds; 
(hi) A221A223, mediating b — > ccd and b — ► sscf; (iv) A / 122 A / 1 23 and A 222 A 223 , giving rise to b — > ccs 
and 6 — > sss (and hence -B^ — > J/ijjKs and 5^ — > (frKs). However, these processes are slepton 
mediated whereas the leptonic and semileptonic B decays discussed here are squark mediated, and 
hence the numbers quoted here should be scaled by (rriq/m^) 2 . The conclusion is that one should 
look for any unusual change in CP asymmetry and/or branching fractions in the abovementioned 
channels to get a supporting evidence for RPV SUSY (though, their absence does not rule out the 
model). 





Ali et al\2t\ 


Greub et al [27] 


Melikhov et al \2£] 


BaBar 


CC: 2.4 x 10~ 5 
CD: 1.7 x 10~ 5 
Inc: 9.7 x 10~ 5 


CC: 9.5 x ICR 
CD: Not possible 
Inc: 2.1 x 10" 4 


CC: 6.4 x 10^* 
CD: Not possible 
Inc: 1.7 x 10" 4 


Belle 


CC: 1.3 x 10~ 4 
CD: 1.2 x 10~ 5 
Inc: 2.6 x 10" 4 


CC: 2.4 x 10~ 4 
CD: Not possible 
Inc: 3.9 x 10" 4 


CC: 2.0 x 10~ 4 
CD: Not possible 
Inc: 3.3 x 10" 4 



Table 6: Bounds on | --^2*2^2*3 1 fr° m the decay B — * Kfi + fi~ with different schemes for evaluating 
the SM branching ratio and two sets of experimental values from BaBar and Belle. CC and CD 
stand for complete constructive and destructive interference between SM and RPV amplitudes 
respectively, and Inc means incoherent amplitude sum. For two theoretical schemes, a complete 
destructive interference is ruled out from the data. 



5. Summary and Conclusions 

In this paper we found the bounds on AA' and A'A' type product couplings coming from 
leptonic and semileptonic r and B decays. Most of our bounds are robust and one to two orders 
of magnitude better than the existing bounds. Some of the updated bounds also have major 
improvement. 

The bounds are only on the magnitude of the product couplings. The phase is irrelevant apart 
from the semileptonic penguin decays of B& mesons. For the latter the most conservative bounds 
come from incoherent amplitude summation but depend sensitively on the theoretical prediction 
for SM BRs. 

After this work was communicated a paper came to the archive where some of these bounds 
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have been discussed ]2J| . 
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